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IN THE SPECIFICATION 

Please replace the paragraph beginning on page 2, line 10 and ending on page 2, line 
13 with the following amended paragraph. 

A positive floating point numb e r vari able x can be represented by an expression 
written as: 

\ x = mx2 e (1) 

where w(l<m<2)isa mantissa and e is a binary exponent. 



Please replace the paragraph beginning on page 3, line 8 and ending on page 3, line 17 
with the following amended paragraph. 

There is therefore provided, in one embodiment of the present invention, a method for 
computing a natural logarithm function that includes steps of: partitioning a mantissa region 
between 1 and 2 into N equally spaced sub-regions; precomputing centerpoints a, of each of 
the N equally spaced sub-regions, where i = 0,. . ., TV -1 ; selecting N sufficiently large so that, 
for each sub-region, a first degree polynomial in m computes log(m) to within a preselected 
degree of accuracy for any m within the sub-region, where m is a mantissa of a binary 
floating point representation of a numben variable; and computing a value of log(x) for a 
binary floating point representation of a particular numb e r x stored in a memory of a 
computing device utilizing the first degree polynomial in m. 

Please add the following paragraph after page 3 and before page 4, line 1 . 




Figure 3 is a flowchart of an embodiment of a method for fast natural log(x) 
calculation. 
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Please replace the paragraph beginning on page 4, line 1 and ending on page 4, line 3 
with the following amended paragraph. 




Figure 3 Figure 4 is a representation of a numb e r variable stored in IEEE single- 
precision binary floating point format, partitioned as in one embodiment of the invention. 



Please add the following paragraph after page 5, line 4 and before page 5, line 5. 



Figure 3 is a flowchart of an embodiment of a method for fast natural log(x) 
calculation. When executing the method, which is described in detail below, computer 36 
partitions 62 a mantissa region between 1 and 2 into N equally spaced sub-regions and 
precomputes 64 a reference point a ( of each of the N equally spaced sub-regions, where 
i = 0,. . ., AT -1 . Computer 36 selects 62 N sufficiently large so that, for each sub-region, a 
first degree polynomial in m computes log(m) to within a preselected degree of accuracy for 
any m within the sub-region, where m is a binary mantissa of a binary floating point 
representation of a variable x. Computer 36 calculates 66 a value of log(x) for a binary 
floating point representation of x stored in mass storage device 38 of computer 36 utilizing 
the first degree polynomial in m , where log(x) is a function of a distance between a, and the 
mantissa. Image reconstructor 34 generates 68 an image by using the computed value of 
log(x). 




Please replace the paragraph beginning on page 5, line 14 and ending on page 6, line 4 
with the following amended paragraph. 

Because (m-a)<l, there are two ways to minimize the error. One way is to increase 
the order of the approximation, and the other is to minimize the distance from m to a. 
Because mantissa m is between 1 and 2, in one embodiment of the present invention, the 
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region between 1 and 2 is partitioned into N equally spaced sub-regions. Centers of each of 
the sub-regions are precomputed and used as reference points in equations (4a) and (4b). By 
partitioning into a sufficiently large number of sub-regions, a low order polynomial function 
produces sufficient accuracy for CT imaging purposes. In particular, by selecting a 
sufficiently large number of sub-regions, for any m within any particular sub-region, log(/w) 
is computed by a first-degree polynomial to within a preselected degree of accuracy within 
that sub-region. For example, computer 36 uses the first degree polynomial in m to compute 
values of log(x) for binary floating point representations of particular numb e rs x stored in its 
memory. 



Please replace the paragraph beginning on page 6, line 9 and ending on page 6, line 26 
with the following amended paragraph. 

Rather than compute a sub-region index using / = round {{m -l)xN), which would 
require six operations, one embodiment of the present invention reduces computation load as 
follows. A partitioning algorithm divides the mantissa of a binary floating point 
numb e r variable in memory into two sub-regions. The sub-regions have index i and Ax, 
where Ax is a distance from mantissa m to reference point a,. Indices / and Ax are directly 
extracted from an IEEE floating-point numb e r variable stored in a computer system, thereby 
reducing computation time and improving accuracy. In one embodiment, mantissa 
partitioning occurs as illustrated in Figur e 3, Figure 4, in which index i ranges from 0 to 127 
and each region represents information extracted from the datum shown in Figure 3. Figure 4. 
More particularly, in a single precision IEEE floating point numb e r, variable, b$i represents a 
sign bit, b 3 o the most significant bit of exponent e, b23 the least significant bit of exponent e, 
bn the most significant bit of mantissa m, and bo the least significant bit of mantissa m, (If it 
is desired to use a different designation for the numbering of bits b, those skilled in the art 
can make the appropriate changes required in the description for notational consistency.) In 
this single precision embodiment, exponent e is extracted directly from bits b$o to b2i\ region / 
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is extracted directly from bits b22 to bj6\ and Ax (a distance from mantissa m to reference 
point aj) is extracted directly from bits bis to bp. 

Please replace the paragraph beginning on page 7, line 1 and ending on page 7, 
line 3 with the following amended paragraph. 

Using the extraction illustrated in Figure 3, Figure 4, a maximum error of 
equation (6) in each sub-region is estimated by an expression written as: 





(7a) 



